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SYNOPSIS 

This thesis deals with dynamic stability of 
multi-span pipes conveying fluid where the flow velocity 
is either constant or a small harmonic component is super- 
imposed on it. The very general equation of motion is used 
to study the stability of the pipes. The element matrices 
are obtained by the finite element (Galerkin) method. For 
steady flow stiffness, damping and mass matrices are 
obtained by finite element m.ethods and a standard dynamical 
matrix is obtained. The critical velocities are obtained 
by solving eigen value problem. For harmonically perturbed 
flow the bounds for the principal primary regions of ins- 
tability are determined by Bolotin's method. It is shown 
that critical velocities and instability regions can be 
controlled by changing the position of the supports. 

Two general computer programs have been written 
to solve these problems. These programs are very flexible 
in nature and solution for various types of multi-span 
pipes can very easily be obtained by changing the input 


data 



CHAPTER - I 

INTRODUCTION 


1.1 INTRODUCTION 

Pipes conveying fluids are encountered in various 
fields of engineering. To name some of them ,* oil pipe 
lines / propellent lines, pump discharge lines, heat 
exchanger tubes and coolant channels of nuclear reactors. 
An imperfect design of pipe lines may result in a leakage 
and unsatisfactory performance of the whole system. As 
the flow velocity is increased a certain velocity is 
reached when the pipe either buckles or flutters. With 
harmonic pertxirbations these systems are vulnerable to 
parametric instability also. 

In the recent years many pipe lines have been 
laid and this has accelerated the research work in this 
field. Lot of research has been done in the last thirty 
years . However most of literature is devoted to either 
single span pipe or periodic pipes. In actual practice 
however we come across nonperiodic multi-span pipes. 

The aim of the present work is to study the 


stability of multi^span pipes conveying fluid at a 
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constant velocity or when a small harmonic component is 
superimposed on it, by finite element method. All types 
of boundary conditions can easily be accommodated by this 
method, .The objective of present thesis is explained in 
detail in the last section of this chapter. 

2.2 REVIEW OF PREVIOUS WORK 

A detailed review of the literature with exten- 
sive references is given in Paidoussis and Issid [21] and 
Chen [sj. A brief review is given below. 

Interest in the subject of dynamics of pipes con- 
veying fluid was developed in 1950 in connection with study 
of vibration of Trans Arabian pipe by Ashley and Haviland 
[1]. They studied the problem of vibration of simply 
supported pipe. The same problem was studied by Housner 
[12] independently by usina a different approach. He 
found that at low fluid velocities the effect upon vibra- 
tions of the pipe was negligible but at certain high velo- 
cities pipe buchled like a column subjected to axial loading. 

Long [15] considered the case of cantilever pipe 
conveying fluid , He adopted an iterative procedure using 
a power series for mode shape, which was applicable to 
relatively small flow velocities and confirmed experi- 
mentally that forced motion of cantilever pipes was damped 
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by internal flow in the range of flow velocities considered. 

Benjamin [2] dealing with the general dynamical 
problem of cantilevered system of articulated pipes con- 
veying fluid was first to anticipate the phenomenon of 
instable oscillations (flutter) . He produced complete 
theory supported by experiments for articulated pipe sys- 
tems . He showed that when the system was vertical both 
oscillation and buckling instabilities are possible. 

However in general when motion is confined to horizontal 
plane / that is gravity is insignificant, buckling can not 
occur. Later Gregory and Paidoussis [lo] confirmed both 
theoretically and experimentally that at sufficiently high 
velocity cantilever pipes are subjected to oscillatory 
instabilities. Next Paidoussis [iP] included the effect 
of gravity and showed that in case of vertical cantilevers 
buckling instability is not possible. This contradiction 
was cleared by Paidoussis and Deskins [ 2 o].. 

Naguleswaran and William [l^] studied 'both theo- 
retically and experimentally the effect of fluid pressure 
on the dynamics of the pipe. They showed that pipes with 
both ends supported may buckle at small flow velocities 
by the action of high internal pressure. 

Chen [5] studied the stability of pipes conveying 
fluid with upstream end fixed and downstream end supported 
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by displacement springs, so that boundary conditions are 
intermediate between clamped free and c lamped -pi nned . He 
showed that both buckling and oscillary instabilities are 
possible depending upon the spring constant. 

Paidoussis and Issid [ 2 l] derived the general 
equation of motion. They included the effect of axial 
movement of the pipe. They studied the dynamic stability 
of the pipe where the flow velocity is either constant 
or a small harmonic component is superimposed on it. They 
showed that with both ends fixed pipe is subjected not 
only to buckling instability but flutter instability is 
also possible at high velocities. 

Singh and Mallik [22] studied the wave propoga- 
tion and vibration response of a periodically supported 
pipe by propogation constant approach. A succeeding paper 
by same authors [23] was devoted to parametric instabi- 
lities . 

Orris and Petyt [is] applied the finite element 
methods to propogation constant in periodic structures . 

Deb [7] applied finite element (Galerkin) method and 
obtained flutter and buckling instabilities of various 
types of single span pipes. 
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Kulkarni [is] studied instabilities of the periodic 
pipes due to parametric excitation considering the effect of 
axial movement of the pipe. He used the propogation cons- 
tant technique . 

1,3 OBJECTIVE AND SCOPE OF PRESENT WORK 

The aim of the present work is to develop finite 
element method for studying the dynamic stability of multi- 
span pipes conveying fluid. 

In second chapter element matrices are obtained 
by the finite element (Galerkin method ) for pipes convey- 
ing fluid when the flow velocity is either constant or a 
small harmonic component superimposed on it. For steady 
flow stiffness, damping and mass matrices are obtained by- 
the finite element method and a standard dynamical matrix 
obtained. Critical velocities are obtained by solving the 
eigen-value problem. For harmonically perturbed flow, 
bounds of principal primary instability regions are obtained 
by using Bolotin's method. 

In third chapter results for various types of 
pipe configurations are given and discussfed. 

Conclusions are reported in the last chapter. 

The values or various matrices and listing of computer 
program are given in Appendixes . 



CHAPTER - II 


GENERAL FORMULATION 


This chapter deals with the finite element 
analysis of stability of pipes conveying fluid at cons 
tant velocity V and with a small harmonic component 
superimposed on it. 

2.1 EQUATION OF MOTION 


Consider a pipe conveying fluid at velocity V 
Fig. 1. Its equation of motion is, Paidoussis and 
Issid [ 21 ] 

E^I _1W — + El -LJi + I M - T + PA (1 - 2vm) 

3X 5 t 8X ■ ^ 


where 


E 

El 

M,, 


- [ (M^ + mp)g - M^ (L - X) [•• 

2 

H- (M^ + m ) -% = 0 



3W 


( 2 . 1 ) 


Coefficient of internal dissipation 
Flexural rigidity of pipe 
Mass of fluid per unit length 
Velocity of the fluid 


V 
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T Longitudinal tension 

P Fluid pressure above atmosphere 

V Poisson's ratio 

^p Mass of pipe per unit length 

L Total length of the pipe 

c Coefficient of viscous damping 

For no axial constraint at the support 
11=1 For axial constraint at the support 

•••• ( 2 « 2 ) 

The equation (2.1) can be expressed in dimen- 
sionless form by defining the following quantities': 


X 

L 


W 

L 


1/2 


+ m 
f P 


1 . 

. 2 




M, 


1/2 M. 1/2 P 

ETMj+nT}^ ^ ' ^f ^p 


7 = 


M_£. + m _ _ 2 

f P T 3 TL 


cL 


PAL 


El 

2 


El 


[EI(M^+mp)]^^^ 


(2.3) 


Substituting equation (2.3) into equation (2.1), the 
equation of motion becomes 


a 


a5 
d w 

4 

3X S T 


d 


+ *1 '-1 - p + p (1 — 2vjij,) + - “K )x 


2 

x(l- x) } u 

3x 


2 2 
^ + r IE + -F Lw 

3x9 T 9X 9ir - 

d V 


(2.4) 


0 
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2.2 PIPES WITH STEADY PLOW 

For constant velocity, equation (2.4) becomes 


a 


— “P + P (l-2vju) - r 

9X 3x 9X gx 


+ 23 


1/2 ^ ^ ^ pi + f «£ 

3x9t >x 9t 


9^7 

9t2 


0 


(2.5 ) 


2.211 Finite Element Analysis: 

For finite element solution (Galerkin Method) 
of this one dimensional partial differential equation, 
the pipe is divided into number of finite elements. 
Figure 2, Typical finite element is shown in Figure 3. 
Its solution is assumed to be of the form of 



or w^®^ ( 5 ^t) = J { ^ (2.6)- 


where 

Interpolating functions of the element 

w. Nodal parameters of the element 
1 / 

r Degrees of freedom of the element (to be 

"decided later) 

5 Local co-ordinate of the element. 


Substituting equation (2.6) into equation (2.5), one 
gets the residue for the element. 
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R 


(e ) 


= a 


jjS (e) .4 (e' 

^ W d W 


W 

4 

rie 9T 


9 w 


+ / -P+ p (l~2vjj.)-r j» 




a2 (e) . 

+ r (X. + 5) + r 


(e) 


+ f 


3 w 


(e ) 


dr 


“f 


H 
2 

9T^ 




1/2 a^w ^ 
+ u ^ ^ 


35 3 T 

(2.7) 


Minimizing this residue by the Galerkin method, 
d5 =0 


or 


g N ( a 


S^w(e) , 


3 5 aor 


35 


4 


2 

+ -^u - n + p (i-2vju) 

(e ) 


. }l5i^ ^ r (X, c) * r Uw 


35 

+ 23^'^^ 


u 




+ f 


3w 


35^ 35 


) de -o (2.8) 


353t BT' Bt;'^ 

In order to reduce the requirements of interpo- 
lating polynomials, the first four terms of equation 
(2.8) are integrated by parts , and one gets 


N. a ( 
1 35^ 


3w^®^ , j 1 

3t 


t 3N. 32 - (e) h h 3 ''n. 

.) I - ^ a I + ; ; 

o 35^ ^*^ 0 0 35^ 


1 33 (e) 

a d^ 

3^i-,3T 


3^6 h 

+ N. — — —■ 

1 353 J, 


3N. 32 ,, (e) h , 3 N, .2 , s 

1 3w i.,h a3w(e)-,.. 

• 1 + / — ^5 d5 


35 


35 


+ .[u - P + p (l-2vM)-r } 


o ~ 35 “ 

(e ) h 


35 


3w 


- {u^ -p + p (l-2vja) - r>-g 


35 ^ 

h ^ 

35 35 


d5 


3w 

+ r X. N. 

1 1 35 


^ V 4. -y M 

i - ^ X,. / q-;: 7 — d? + N 


•35 


35 


(e ) h 
i 3~ o 
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’k 5N. - (e) , 

_ T, fh a ?!w ^ .h 

^ o a? ? 3E i N. 2P 

X 


+ f K 3" 


B? ^ as 

(e ) 


V2 ^ iVl! 35 
55 3T^ 


3t 


Vi p,2 (e) 

d5 + ^ N, — dC = 0 


31 :' 


(2.9) 


(e ) 

Highest derivative of w with respect to g in any 
integral in equation (2.9) is second. Thus interpolat- 
ing function N. should have compatibility of and 


3w 


(e) 


35 


Highest derivative in equation (2.9) is third. 


Thus interpolating function N. should have completeness 

(^\ a,,(e) . 2 ,, .3 (e) 

, ve; 3w 3 w ^ d w ^ 

ox vf , — -- — , — ^ and = 5 — . The function satisfy- 

95 35-^ 353 

ing these requirements of compatibj lity and completeness 
is 

( 0 ) 2 

w = a + bx + cx + dx"' 


or w'®' = i 5 ^ ( 3 - 2 Cj),hCj 52 


4 0-2 2),-hegh 



or w 


(e) 


j_N (5)J { W ( r )J. 


ne 


where 5^ = 1 - 5 A/ 

that is the typical element has four degrees of freedom, 


Using eauation (2.10), the equation (2.9) can be 


written in the following matrix form 
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1 N a 

^ ) j 1 


a 

^2 

(3w^ V 

1. 

oi 

‘ 1 


3^2 

i 




33 , ( e ) h 


w 


,J 


H' 




r 2 

< N(u -f + p (l-2v]U) - 7 ) 

i 


8w 


o 

(e) h 


° i 


-fy N (x. + 5 ) — , 

J 9C J 


(e ) ■] h 
o 


( 2 . 11 ) 


where 
(e) 


[m] 


? {n} L n j as 


(2.12) 


= a N" JdC + 2i3^/2 ^ /^{NHN'JdS 

^ ^5^ {^ } f ^ i (2.13) 

N^'JdC 

- (u^-p + p(l-2''ju) - r) { N' J[n' J CU 


“ {^'} J d5 - r ^5 ^{n'|^ N'j d? (2.14) 

Values of matrices and [k]*" ^-are 

given in Appendix I. Using the values of Ih from equa~ 
tion ( 2 . 10 ) in the right side of equation (2.11), one 


gets 
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+ [c:](e)_^^yne) ^ j-^j (e ) ne ) 


4*4 


a — L 

a r 3 T 

"1 

I 


0 


r 3 1 

3 w ' j 

h- 




35 ^ f j 

! 0 



„ 3 ^ 7 ^®^ } 

- a 


0 i 

f ^ 

■ 3^3 aT 

1 . 

> + < 

axj 

0 

< 

1 + 

1 

[ 

! 

A<e)j ;' 
3 5 ^ ^k 

0 



^3 ( e ) : ; 

„ 9 w ! 

Ct. 2 •: 


1 

0 1 


J 


L. 3 5 3 ^ k J 


1 

Jf 


n 


O 


.. I. 


0 




(u^- r +p(i-2vu)-y) 


0 


B w 


(e) 


3€ 


►+‘)-(u^-r +p(i-2vw)-r) — I 
I • 9 5 lx 


rx- 


5 w(e ) 
3C 

0 


■% 




(e) 


3r 


0 


35 


0 


• • • ( 2 « 15 } 

The equation (2.15.) is for one finite element. These are 
assembled for the whole domain, 

[h] ( w + [C] { + [k] j wp") = 
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i ^ 3\/e) 

i 

r 

0 

t 

Ak<®>i 

35^ 



3 

Br 1 

0 

J : 

; 

1 

' 

+] 

^3 f e ) I 
- a ^ r 

35^ 3T 

# 

^ i 

'"-f < 

0 

• 

» 

-a » 

9 5"* 3“^ !n+l 


m ) 

0 


^3 . 

_ __J1 

3 5^ ' n+1 

0 


CY _ j * 

. 2 t ^ 


0 

- i 

i 


0 




Hr 


0 


{u -y +p(l-2vM)~r ) 


aw 


(e ) 


3? 


O 


.2 •. 
9 W 


ar 


+<- (u^-iT +p(l-2v/x)-r ) 


(e )i 


+< 


■y rih 


n+1 

0 

O 


3w^® ^ 


35 


ii+1 


0 


35 


n+1 


(2.16) 


,0 


After applying boiindary conditions, the equation (2,16) 
reduces to 
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I (n) 


[c] 


BC 




1 (n) 


[k] 


BC 


/ ■ 
] rf 


(n) 


For boundary conditions, see Appendix 2. 
2.2.2 Method of Solution; 


0 

(2.17) 


Equations (^^.17) are set of hornopenous differen- 
tial equations . The solution of this homogenous set of 
differential equations is obtained as detailed in 
Meirovitch (I 4 ]. But it may be noted that the matrices 

symmetric as stated 

in Meirovitch [I 4 ], see Frazer, Duncan and Collar [ 9 ]. 
Using generalised velocity vj as auxilic.ry variables, 
n second order ordinary differential equations (2.17) 


are converted to a set of 2n first order ordinary diffe- 


rential equations, 

= jo] 


where 


^ ^ (%) 

[ . i 


and 


[m] = 


■L 


W (T ) 


[ [ " 


[ 0 ] [m] 3 , 

^^^BC^'“^BC 


and 


[k] = 


[o] 


(2.18) 


(2.19) 


[o] 


( 2 . 20 ) 
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Now we have the differential equations which in 
the matrix form are 


[m] + 

[k] { 

n(T) } 

II 

o 

(2.21) 

set = 

w 

e" 

1 ^ F 


(2.22) 

where 0 consists of a 

. vector 

consisting of 

2n constant 

terms . Substitution of equation 

(2.22) in equation (2.21) 

leads to the eigen 

-value 

problem 



1 1 

j 1 

1 1 

+ 

1 } 

= 0 

(2.23) 

which can be written as 




[d] ( 0 'f = 

1 

f 1 

0.> 


(2.24) 

where 






1 

1 

[0] 

[l] 


(d] = - [k]-^ 

It 

t i 

1 1 

1 

I 

m] -[k]"^[c] 

(2.25) 


The eigen value of equation (2.23) will in general 
be complex. The real part of eigen value determines the 
stability of pipe and imaginary part gives the frequency 
of pipe. If the real part is positive and imaginary part 
is nonzero then there will be flutter- If real part is 
positive and imaginary part is zero the pipe will buckle. 
The velocity u at which pipe fails either due to buckl- 
ing or flutter is known as the critical velocity, u^. 

A library subroutine EIGRP was used to calculate 
the complex eigen values . 



2 . 3 PIPES WITH 


HARMONICALLY PERTURBED FLOW 


Equation (2.4) is general governing differential 
equation of the pipe conveying fluid. 


Now a case is considered where the flow velocity 
is perturbed harmonically, that is 

u = u^ (1 + 6 cos Qt ) (2.26) 

where u^ is mean flow velocity, 6 is excitation para- 
meter and 3 is non-dimensional frequency given by 

*ir 

Substituting equation (2.26) in equation (2.4) 

one gets 

5 4 

O', — - . 3h , + ^ Qt )^ P (l-2v|i) 

rtx Bt BX 


-r (1 - x) - 6q (1 - x) sin Q t] 

9x 


20 1/2 ( 1+6 cos Q T ) 

o 

,2 

+ = 0 


' n 2 „ , 

9 W . .y ^ 

Bx BT ^ Bx ^ BT 


BT^ 


(2.27) 


The coefficient of some derivatives are time 
varying, hence the pipe is parametrically excited. Because 
of parametric excitation, instability will occur over a 
range of frequencies. 


In order to obtain the primary instability 
regions, we can write, Bolotin [ 4 ] 
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w ( X, ) = S (Y^ (x) sin (^) + (x) cos (^) 

k = 1, 3,5 (2.28) 

The regions of primary instabilities can be 
obtained quite accurately by taking K = 1 approximation, 
Paidoussis and Issid [2l]. Hence one can write 

w ( X , tt ) = Y (x) sin + Z (x) cos ( - ■ — ■ » ) (2.29) 

Substituting equation (2.29) in equation (2.27) 
and equating the coefficients of sin (^^^-^) and cos ^ 

separately to zero, one gets 


^ (1-6+ _p + p (1 _ 2viu) - r (1-x))^^ 


dx 


dx 


R u 6 p 2 

1 ^2 _ dY ^ ^o r/-, ^ d^Z 2, dZi 

a y + r gj 5 [(i-x) ^ - (1 - J.) g-J 

dx 


a Q d Z f S 

2 , 4 " 2 
dx 


0 


(2.30a) 

.2 


4 2 

^ ^ -f ( u (1 -f 6 -f. ).«ji ‘fp(l — 2 v/x ) — 7 (1 --x)) 

dx dx 


6Q 2 

I ^ " i - ^ ! — hi - X) ^ - (1 + =) g] 

dx 


, n Q dS" . f Q 
+ ^ + - 5 - Y = 0 

dx 


5-' d^- 
( 2 . 30b ) 


Let us take 


A 

A 


1 

2 




< 5 +^)-p + p(l-2v )-r 


(2.31) 


(2.32) 


MINI 
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u 




^4 = 


OQ 


(2.33) 

(2.34) 


^s-u^ (!+£+—) - +p(l - 2viJ.) - r (2.35) 

8^”^^ U 6q 

Ag = (1 + -) 2 2 (2.36) 


Using equations (2.31) to (2.36), equations (2. 


become 


,2. 




d^Z 

dx" 


+ A 


dZ 

a Q 

dS _ 

f Q r. 

^4 ^ 

2 

d^Z 

dx^ 

— 2 

d 

; ^5 

5? " 

A 2 2 

+ -y ^ 
dx 

dY 

. a Q 

d'^Y 

f Q ,, 

6 dx 

+ - 5 - 

dx 

— Y 


0 


/ ( 3.2 % 

(x 3-) 


A3 (1 


- X) 


(2.37a) 

d^ 


dx' 


(2.37b) 


2.3.1 Finite Element Analysis: 

Equations (2.37) are coupled differential equa- 
tions. For finite element solution of these two coupled 
differential equations, the pipe is divided into number 
of finite elements. Figure 2. Solution over the finite 
element is assumed of the form of 


^(e) 


! N . ( s ) i J Y . ^ 

L ' ^ ..i L 1 J 


30) 


(C) = 


(2.38a) 
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and 

( 5 ) = IB. (?) i j Z (2,38b) 

. 1 j 

«r ^ 

where , 

Interpolating functions 
Nodal parameters 
5 Local coordinate 

Substituting equation (2.38) in (2.37) one gets the 
residues ^ and for equation (2.37a) and (2.37b) 

respectively. 


( e ) d ^ ^ ^ (3 ^ Y ) a 

+ A A. + r ^ ((x. t?) -^1 ) 

^ d?"' ^ d?2 2 d? J d? 


A, (1 - X, -5) + A, 22 A - ^ 2 '®' 

3 J ■ 4 d; 2 ^^4 2 


d^z^"^ . , dz 


( 2 . 39 a) 

R^e) ^ + A + r ((x +?) '^1^) 

d£ ^ dP^ ^ d? J d 5 


. A, (1 - X. -c) + A, + 22 fO ^(e) 

3 J ds^ ® -35 2 3^4 2 

( 2 . 39b ) 

where x ^ is the global coordinate of the j node . 

Minimizing -the residues by the Galerkin method 


K. R^'®> as = 0 


(2.40a) 




o 


(2.40b) 
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Treatment for (2.40a) and (2.40b) is similar. 

So details are given only for (2.40a) and results of both 
are given in equation (2.46). 


(f=‘) 

Substituting the value of R “ from equation 

3l 

(2.39a) in the equation (2,40a)/ one gets 


N. (^-I^ — + A. A, + r ((x.+ 5)~ 

6 1 dS^ 1 ^^2 2 dC 1 dC 


) 


, , dZ^®^ a Q dS^®^ f Q 

- ^3 + ^4 ^ T- 774 ^ ^ 


dC 


dS 


0 


2 

(2.41 ) 


In order to reduce the requirements of the enterpolating 

polynomials first, second, fourth, fifth and seventh terms 

are integrated by parts, and one gets 

, d^N. ,2 (e) T_ dN. n„(e) , 

i ar^ as - Aj g'’ gU «_ d6 - N^Y ai 

- g c! ar if- '"s - M E 31^ i|- as 


+ A (1 - X ) I 


h aN^ a^Ce) 


, dN . dz 

h 1 


(e) 


j' o d^ d^ 
(e) 


“ ^3 ^ d? d^ 
(e) 


d5 


.h „ dZ'"''. . , ,h dZ 

- A^ / N.. :3-r — d^ + A_^ / d^ 


’3 6 1 d? 

, d^N. ,2 (e) 
aQ ,h idz 

i - 3 - — 7 - as 

d 5 d 5 
3„(e) dN, 


^ N. Z^^^ d? 

26 i ^ 


N 


d-2'-^ h d^^®^ h , , dY^®^ h 

= - — - — K 1 + A- N. ■3T ' 

^^72 ! 1 1 d5 


^ dr 


i 


dC 


6 
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- A3 (1 

- X . - 

1 

5 ) ■ 

i 

- T B. (Xj + 

,, dY'®’ 

0 

ry, Q 

ah<«' 

dC^ 

h 

1 + 
6 

a B 

2 dS 

dS^ 0 

(2.42a) 


Similar equation (2.42b) is obtained from equa- 
tion (2.40b). Highest order of derivative in the inte- 
grals in equations (2.42) is two. Therefore and 

( p ) ^ Y V 

Z ^ should have compatibility of (Y, ^) and (Z, g^) 
respectively. Highest order of derivatives in equations 
(2.42) is three. Therefore y^^^ and should have ^ 

. ,2 Y 

completeness of (Y, ^ , — ~) and (z, ^ 

d?^ dC^ dY dC^ 

respectively. Thus and can be chosen identically. 

Interpolating functions satisfying these requirements 

are. Able and Desai [s] 


|Nj_(?)j»jB. (5)t= \jl (3 - 2z),h^lK2 ^ \ 

v/here ~ ^ ~ h ^2 ~ K (2.43) 

Using equation (2.43). the equation (2.42a) 


becomes 


- *2 }U J 

- •' tN’ J 


h } )Y'J “SS 

- ■>' a5{Y}'“> 

+ Ajd - Xj) 
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^3 J' • J /^{N '>jN' i dg {zl 


+ ^4 '! N } |N' Jdc '[“zl^^®^- ^ N'' I In'' jdC <i_z| 


I ( ne 





: 

"■ 1 

i ^p3 

d C J 


0 

0 


d^Y^®^ 1 

— ^ ' ' ■ 1 

(i 


de^ , 

< 

+ < 

1 

d\(e)v 1 


f 

^ i 

d?^ 'k 


0 

0 


! dk‘«>, 

u -• 

i- ^c2 T, J 


dY 

3J 


(e) 


0 


+ Ai » dY^®,- 




k 



~- (1 

n ( 0 ), 

-"’dl } 

1 

J 



0 


+ A- •«; 
I" 3 

(1 

0 ) 

N W' 

1 

. 

' 

k 

J 

i 

\ 

C. 

0 



r (e ) 

dY^ I 

-- X 


+ r 


dg 


O 


X 


dY 


,(e) 


dK 


0 


>+ <■ 


k 


d^Z^®! 
dr^ ‘ 


0 


'I 

dX~ 


0 


0 




ds' 


+ £_£ < 
2 


0 


I dS^ 


k ^ 


(2.44a) 
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Similarly equation (2 •42b) can be written as 


t { f’’ 


dc (z) 


(ne) 




dC / Z 


(ne ) 


i- J 


- A, /“ J n1 

2 O L J 


N|dc r X d{N'>lN'jd5 {Z}'"®’ 


J o i 


>lN'Jd5 {z} 


y 6 5'[ N'} L N', j dC Z ■ 


(ne ) 


+ A^(l-x . ) i N'l I N' (d? 


(ne ) 


- Aj dsl N'} L H'J dE J Y 


■3 ^ {N rlN'J dE JY} 


i ( ne ) 


+ Ag { N T,p' J d5 ^ J d£J 


+ ^ r -iN-||N-Jd€|Y}' 


V L n j d5 [y^ 


1 ( ne ) 
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25 


+ 


O’. Q 






J '> 


O 




dF/ 


i 


0 


d^Y^® ) 


J 

N 


d5‘ 


0 


+ 


d\^® ) i 
dl 


kJ 


(2.45) 


Matrices [a], [b], [e] and [g] are given in 
Appendix 3 . 


The equations (2.45) are for one element. Now 
assembling it for whole domain and applying boundary 
conditions , one gets 


or 


’[a] 

[b] 



' {- 7 1 

[e] 



{ 1 


^ n 



i 

I — 1 

= 1Y/ 

Jr = 


0 







Various types of boundary conditions are listed 


in Appendix 4 
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2.3.2 Method of Solution: 

Bolotin has shown that instability regions can 
be obtained by equating determinant of matrix [p] to zero 
Hence frequencies Qj , j = 1, 2.,.,. are the frequencies 

where Det [f] :=; o and they give boundaries of instability 
regions . 

A Library subroutine F03AAF/NAG was used to cal- 
culate the determinant of matrix [f]. The instabilities 
regions are plotted in Q 6 space* 



CHAPTER - III 


RESULTS AND DISCUSSION 

This chapter deals with numerical results obtained 
for pipe conveying fluid for various types of pipe confi- 
gurations. In Section 3.1 stabilities of pipes are dis- 
cussed when the fluid is flowing at constant velocity. 

In Section 3.2 instability regions for various types of 
pipe configurations are plotted when a small harmonic 
component is superimposed over the mean velocity. 

3.1 PIPES WITH STEADY FLOW 

In this section critical velocities of various 
types of pipe config^Irations have been found by the finite 
element methods. 

3.1.1 N\imber of Finite Element: 

To begin with one must decide -tiie number of finite 
elements to be used in the analysis. For every type of 
boundary condition the frequencies were calculated for 
few velocities with increasing number of finite elements. 

We selected that number of elements beyond which the 
increase in number of finite element did not change the 
value of frequency significantly. 
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Table 1 shows the frequencies of two equal span 
simply supported pipe with increasing number of finite 
elements for =0.5, u = 2 and a=;'=p=r = f = 0. It 
is observed that the value does not change significantly 
after eight elements. So the number of finite elements 
taken for this case was eight. 

3.1.2 Results and Discussions: 

Using the finite element equations developed in 
Chapter-II the critical velocities for various configura- 
tion of pipes are obtained. 

3. 1.2.1 One Span Pipes 

Complex frequencies for one span simply supported 
pipes were obtained by Paidoussis and Issid [2l] for 
3 = 0.5, a=p’=p = r = f = 0, The same problem was 
studied by Deb [?] by using finite element method. 
However Deb [?] did not study the effect of internal 
dissipation (a) and gravity (r) effects. 

Figure 4 shows effect of coefficient of internal 
dissipation on the same system. The value of a is 
taken equal to 0.005. It is found that critical 
velocity does not change but the symmetry about 
imaginary axis is lost. These results match with 
those obtained by Paidoussis and Issid [2l]. 
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critical velocities of all these four cases are less 
than the critical velocity of three equal span pipe, 

3. 1.2. 3 Multi-span Pipes, One End Pi and Other 
Supports, Simple Supports 

Critical velocity for tv/o equal span pipe, case (1) 
Table 3, for .6 = 0.5, a =p= p = r = f = 0, is 7.12. 
This critical velocity is more than the critical 
velocity of pipe v;hen all the supports are simple 
supports as expected. Next this pipe is considered 
with two unequal span lengths, cases (2, 3), Table 3. 
It is found when the intermediate support is shifted 
away from the fixed end critical velocity becomes 
7.42 which is more than that of the pipe considered 
in case (1). On the other hand if the intermediate 
support is shifted towards the fixed end the critical 
velocity becomes 5.56 which is less than that of 
case (1). 

Next three span pipes are considered, case (4, 

5, 6,7), Table 3. It is found for equal span pipe 
the critical velocity is 9.93. Critical velocities 
for unequal spans, cases (5, 6, 7) are 8.60, 8.10 
and 10.06 respectively. It may be noted that these 
critical velocities are more than those if all the 
supports are simple supports, see Table 2. 
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3. 1.2. 4 Multi-span Pipes with One End Fixed, Other 
End Free and Intermediate Supporcs Simple 
Supports 

Critical velocity for this tv/o equal span pipe, 
case (1), Table 4 is 4.64 for P =0.5, a=P=p = r 
= f = 0. It the position of intermediate support is 
shifted towards the free end the critical velocity 
becomes 5.65, case (2), Table 4^ On the other hand 
if intermediate support is shifted toxvard the fixed 
end, the critical velocity becomes 3.96 which is less 
than that of pipe considered in case (1). 

Next three ■ span pipes are considered, cases 
(4, 5, 6, 7) Table 4. The critical velocities for 
these cases are 6.48, 7.41, 5,79 and 5.24 respectively. 
It may be noted that these pipes are similar to these 
considered in cases (4, 5, 6, 7) of Table 3, except 
that the right hand support has been removed and is 
free now. It is observed that if the right hand 
support is free, the critical velocities decrease. 

3. 1.2. 5 Miilti-span Pipes with Both Ends Fixed and 

Supports Simple Supports 

Critical velocity for two equal span pipe, case 
(1), Table 5, for 3=0. 5, a=p = p=r = f = 0, is 
9.0. For unequal span lengths case (2), Table 5 the 
critical velocity decreases to 8.65. 
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For three span pipes, cases (3, 4, 5, 6) critical 
velocities are 11.76, 9.90, 11,24 and 10.26 respec- 
tively . It should be noted that these pipes are 
similar to those considered in cases (4, 5, 6, 7) 
Table 3, except that the right hand support has been 
fixed. It is observed that in the present cases 
critical velocities increase. 

3. 1.2. 6 Multi-span pipes with One End Fixed 
and Other Supports Elastic Supports 

Table 6 gives the results for the pipes whose one 

end is fixed and all other supports are elastic 

supports. In cases (1, 2, 3, 4), Table 6 elastic 

supports are displacement springs. For 3 = 0.5, 

a=p=p=r = f = 0 and spring constant = 10. o 

the critical velocities are 8.52, 7.14, 7.08 and 

6.93 respectively. These velocities are less than 

that of those considered in cases (4, 5, 6, 7), 

Table 3. Where all the supports are simple supports 

as expected . 

Next pipes with torsional springs are considered, 
see cases (5, 6,7, 8), Table 6. The critical velo- 
cities are 13.41, 10.92, 10,71 and 10.38 respec- 
tively. These velocities are more than that of these 
considered in cases (4, 5, 6, 7), Table 3 as expected. 
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3. 1.2, 7 Multi-span Vertical Pipes 

Table 7 gives the critical velocities for vertical 
multi-span pipes. It is observed that critical velo- 
cities are more in case of downward flow than those 
for horizontal flow, and critical velocities are less 
in case of upward flow. It is observed that effect 
of gravity is more pronounced in multi-span pipes 
if one of the ends is free . 

3.2 PIPES WITH HARMONICALLY PERTURBED FLOVJ 

In this section we will discuss the stability of 
pipes when a small harmonic component is superimposed over 
the mean velocity. 

3.2.1 Results and Discussions: 

Using the finite element equations developed in 
Chapter-II the regions of principal primary instability 
are obtained for various types of pipe conf igxrrations . 

These regions are plotted in Q - 6 space. 

3. 2. 1.1 One Span Pipes 

In Figure 5 instability regions for one span 
simply supported pipe are plotted for u = 0,6t[, 
a =p = p= r = f = 0 and = 0,4 and 3^'^^ = 0.8 

These results match with those obtained in [21,13]. 
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1/2 

For 3 =0.8 results are also given in Table 8 and 

compared with Kulkarni [l3]. Excellent matching is 
observed. It is observed if the value of 3 is increased 
instability region broadens. 

Results for fixed-fixed pipe are also given in 
the Figure 5. Instability regions are plotted for 
3^-^^ = 0.4, a=p=p=r = f = 0, u = 0.6rr and tt. It 
is observed that with increase in velocity u the 
instability region shifts to lower value of Q. 

In Figure 6 instability regions are plotted for 

1/2 

fixed simply supported pipe for 3' =0.4, a= 

= p = 7 ' = f = 0, u = 0.6'’T and rr. Again it was 
observed that with increase in value of u instabi- 
lity region shifts to lower value of Q. In the same 
figure instability regions for cantilever pipe are 
also plotted for 3^'^^ = 0.4, a =P =h= r = f = O, 

■u = m and u=1.2Tr. A striking feature of this 
instability region is that it starts after a finite 
value of 6 , Again instability region shifts itself 
to lower value of Q for higher velocity. 

3. 2. 1.2 Multi-span Simply Supported Pipes 

Instability regions associated with first two 
modes of two equal span simply supported pipe for 
u = 1.2 7T, 3^^^ = 0,8, a=p=p = r = f = 0 are 
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given in Figure 7 and Table 9. The results obtained 
by Kulkarni [is] are also given in the same ■ figxire 
and table. By present analysis I'le obtained broader 
instability regions. The difference arises due to 
the change in the governing differential equation of 
motion. Kulkarni [l3] has taken L in third term in 
equation (2.1) as length of the si in, X'/here actually 
L is the overall length of the pipe. 

Figure 8 and Figure 9 give the instability regions 

associated with first two modes -of tx*;o and three span 

1/2 

(equal and unequal) simply supported pipe for ^ =0.4, 

a=n =p=r=f=0, respectively. For two span 
pipe velocity u is 0.6 r and for three span pipe it 
is 1.8 rr. It is seen that for unequal span lengths 
the instability region associated with first mode 
shifts to lower value of Q in comparison to first 
mode of equal span pipes . On the other hand the 
instability region corresponding to second mode shifts 
to higher value of Q in comparison to instability 
region corresponding to second mode of equal span 
pipes. The instability region corresponding to first 
mode of one span pipe is also plotted in Figure 8. 

It can be seen that for one span pipe the instability 
region are at very low value of C. 
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3. 2. 1.3 Multi-span Pipes, one End Fixed and 
Other supports Simple Supports 

Figure 10 and Figure 11 show the instability 
regions associated ivith first two modes of two and 
three span (equal and unequal) pipe for 3 ^ = 0,4, 
a=P = p=:y = f = o respectively. For txvo span 
pipe u is 1.2 rr and for three span pipe it is 1,8 n. 
It is observed that regions of instability associated 
with first mode of unequal spans are at lower value 
of Q in comparison to that of equal spans. On the 
other hand the instability regions associated with 
second mode are at higher value of Q in case of un- 
equal spans. In Fig\jre 10 instability regions corres- 
ponding to first mode of one span fixed - simply 
supported are also plotted. It is observed that for 
one span pipe instability region is at very low value 
of Q . 


3. 2. 1.4 Multi-span Pipes, Both Ends fixed and 
Other Supports Simple Supports 

FigTxrel2 gives instability regions associated 
with first two modes of a two span (equal and unequal) 
pipe for u = 1.2 tt, = 0.4, a=p = p=r = f = 0. 

The region associated with first mode is at lower 
value of Q for unequal span length. But region of 
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instability associated with second mode is at higher 
value of Q for unequal spans . The results for one 
span are also plotted for comparison purposes . 

Figure 13 gives instability regions associated 
with first two modes of three span (equal and unequal) 
pipe for u = 1.8 n, 8^/^= 0.4, a ='''-=p = r = f = o. 
Here the instability regions associated with both the 
modes are at lower value of Q for unequal spans. 

3. 2. 1.5 Multi-span Pipes, One End Fixed, Other and 

Free and Intermediate Supports Simple Supports 

Figure 14 shows the instability regions for a 
two span (equal and unequal) pipe. The interesting 
feature is that the instability region starts from 
a finite value of 6 and both modes are at higher 
value of Q for unequal spans. 

Figure 15 gives instability regions for a three 
span pipe for u = 1.8 rr, = 0.4, a =n =p=r=f=0. 

The instability regions of both modes are at lower 
value of Q for unequal spans. The instability region 
corresponding to the second mode brodens in case of 
unequal spans . 



CHAPTER IV 


CONCLUSIONS 

Based on the results obtained in Chapter III folloxv- 
ing conclusions are drawn. 

1 . The results by finite element method for multi- 
equispan simply supported pipes with steady flow match with 
those obtained by earlier authors. This shows that finite 
element equations of motion are correct and can be applied 
for multi-span (equal and unequal) pipes with any type of 
supports . 

2. The results by the finite elements method for multi- 
equispan simply supported pipes with harmonically perturbed 
flow match with those obtained by earlier authors for P - 0. 
For P ^ 0 results differ from earlier results because of 
correct equation of motion used here. Thus the finite ele- 
ment equations of motion developed here are correct arKl can 
be applied for multi-span (equal and unequal) pipes with any 
type of supports . 

3. The method developed is very general and flexible. 
The critical velocities and instability regions for any type 
of pipe configurations can be obtained by changing the input 
data in the computer program. 
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4. The column matrices! N(u‘^ - p +p(l-2vM,) - r) ^ 1 f 

and V N(xj | i’ of equation (2.11) should be trans~ 

ferred to square matrix of left hand side, similarly, third, 
fourth fifth column ma tries on right hand side of equation 
(2.45) should be transferred to square matrices of left hand 
sides, as suggested by Deb [?]. These column matrices can 
not be neglected as in suggested by the work of Szabo and 
Lee [ 24 ]. 

5. Critical velocities of pipes can be controlled by 
changing the position of the supports. Critical velocities 
of vertical pipes decrease for upward flow and increase for 
downward flow in comparison to horizontal flow. This 
effect is more pronounced if one end is free . 

6. The principal primary instability regions broadens 
when the value of 3 is increased. The principal primary 
instability regions shifts to lower value of Q for increase 
in velocity u. 

7 . If one end of the pipe is free the , principal 
primary instability regions starts from a finite value of 

6 

8. Instability regions can be controlled by changing 
the position of the supports . 
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TABLE 2 
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TABLE 5 

velocities for Multi-span pipes with boLh end fixed 
and other supx^ort s sim ple supports 

P = OoS, a=p = p = 7 ' = f = 0 
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I 


TABLE 4 

Critical velocities for Multi-span pipes xvith one end fixed / 
other end free and intermediate supports s imple s\ipports 

P = 0.5, a =r* = p= r = f= o 





50 


TAF’-S 6 


Crj-bi .ca 1 velocities for Multiple span pipes v ; Ith one end fixed 

and other supports elastic supports 

,G = 0.5, a = p = p = 7' = f = 0 
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TABLE 7 






TABLE 3 


Values of frequencies bounding principal 
primary instability region for single span 
simply supported pipe 

u = 0.6 n, = 0.8, a =n = p = r = f = 0 

NEL = 5 



Freer 

lencies 

Lo^\rer 

i Upper 

0.01 

15 .345 

15 .507 

0.10 

14.613 

16.204 

0.20 

13.765 

16.948 

0.30 

12.883 

17.654 

0.40 

11.955 

■ t 

18 .319 
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TAELS 9 


Values of f recTTgencies bounding tihe principal primary 
instability regions associated with first two modes 
of a two span pipe obtained by present analysis and 

by KulTcarni flsl 

u = 1.2Tr, = 0.8, a=P=p = rf = 0 

NEL = 8 


Fxrst Mode Second Mode 


Lower 


sent 


ent i karni < sent karni 


0.01 

61.20 

0.10 

56.07 

0,20 

50.38 

0 . 30 

44.87 

0.40 

i 

1 

j 

39.60 


0.01 61.20 61.40 j 62.29 62.40 | 109 . 60; 109 .88 110 .92i 110 .44 

i 

0.10 56.07 58.52 67.17 65.72 104.561 106 .08 116,67 113.06 

0,20 50.38 55,02 71.97 66.56 99.69 104.04 123.84 116.08 

0.30 44.87 51.60 75.59 71.12 96.02 lOl .00 1 3l .50 118 .84 

0.40 39.60 47.88 177.10 72.36 94.21i 98.16 139.241131.48 
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Fig 2 Ffnlt« tiemcnts of tht pipt 



Fig. 3 Typical finite element of the pipe 





U = 0.6TT>P " 


^ = 0‘i 

sivp. — — * 


0.6 %P 





Regions ot principal primary instability of 
one span pipes for 

oLrr-.p- J=:f=o NEL = 5 



0.1 


0.3 


0 . 


0.2 

S 


Fig. 6 Regions of principal primary instability of 
one span pipes for 

0.4> oCrPspsit faO NELa 










Fig. 8 Regions of principal primary instability of a 
two span simply supported pipes for 

tta 0*6 TT^- 0*4 > eta Fa pr”? a faO NEL. a 8 






Fig. 9 Regions of principal primary instabHIly ®f thi 
span simply supported pipes for 

U*1.8Tr, P^^^a0.4>ol.*raps**t»o NEL*12 



6 

FIg.tO Regions of principal primary instability of two 
span pipes one end fixed and other supports 
simple supports for 
U=1,27r,p’''=0.4,ot=r=P=-i-f=o NEL = 8 







simpi® support for 

u=t. 2 Tr, p’'^=0.4.ot=r*p=i = f = o NEL = 8 



Regions of pincipol primary 

SPOT ptp« both ends fixed and other supports 

simple supports wpi -i? 

u=i.8Tr,p’''=o.A,<i:=r=P=#*f=« ne'--i2 



Ffg,U Rcgfons of principal primary instability of two 
span pipes one end fixed other free and 
Intermediate support simple suppot for 
u=1.2TT,p''^:.0./.,o( = r=p=V=f = o 


NEL = 8 




Fig. 15 Regions of principal primary instability of thr^ii 
span pipes one end fixed other free and 
Intermediate supports simple supports for 

UsI.aiT# P^^^s04»<iarspal»f«o NELcl2 
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APPE^TDIx-I 

The values of various matrices are as follows 


r 

i.mj 



156 

22h 

54 

- I3h 

h 1 
420 


4h2 

I3h 

-- 3h^ 

1 



156 

- 22h 

! 


sym 


CM 


(Al .1 ) 


[c] 


(e ) 


or. 


h' 









— 

12 6h 

- 12 

6h ' 


-30 

6h 

30 

-6h 

4h2 

- 6h 

2h^ 


-6h 

0 

6h 

-h2 




^ 2!3l/h 






12 

- 6h 

60 ■ 

-30 

-6h 

30 

6h 

sym 


4h2_ 


6h 

L 

h2 

-6h 

0 

i 


4 - 


f h 
4 ^ 


[k] 


(c) 


1 


12 


156 

22h 

54 

- I3h 


4h2 

I3h 

- 3h 

sym 


156 

- 22h 




4h 

6h 

-12 

6h 1 



(A1.2) 




4h' 


sym 


-6h 2h 

12 -6h 

2 


+p(1^2vji)~-y) 

” 30h 


36 


sym 


4h 


4h 

3h 

2 


36 

3h 

36 


3h s 
2 


- h 


4h 


3h 

2 


Contd - 
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36 

3h 

- 36 

3h 


36 6h 

- 36 

0 

r X . 

L. ' 

30h 


.1^2 

4h 

- 3h 

v2 

~h 


2h^ 

- 6h 

-h^ 



36 

-3h 

60 


36 

0 


sym 



.•u2 

4h 1 


sym 


CM 

CO 





J 






(A1.3) 
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APPENDIX- 2 

BOUNDARY CONDITIONS FOR STEADY FLOW 
A 2.1 Simple Supports : 

For pipe simply supported at extreme ends, the 
boundary conditions will be 


w = 0 


3^2 


0 


(A2.1) 


For a intermediate simple support, the boundary 
con<lltion will be 

w = 0 (A2.2) 

A2.2 Fixed Supports : 

For fixed ends of the pipe, the boundary conditions 

will be 


w = 0 


aw 


= 0 


(A2.3) 


A2,3 Flastic Supports, Displacement Springs ; 


For pipe supported by a displacement spring, 


+ a-w = 0 atC=0 
3 w 


w = 0 at 5 
d 


(A2.4) 


where 


a 


K^L- 


'd B I 


(A2.5) 
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A2.4 Elastic Supports, Torsional Springs : 

For pipe supported by torsional spring, 

9 


3 

t 35 

= 0 at 

ii 

3 w 

'Ip 

+ n. 
t 35 

= 0 at 

\fn 

ii 


where 


a. 


t 


K.L 

t 


V 


(A2.6) 


El 


(A2.7 ) 
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appendix ~ 


The matrices [a], [b], [e] and [f] can be given as 


[a] = 


12 6h -12 


6h 

2 


A, 


420 


4h -6h 2h 
sym 12 -6h 

4h 

156 22h 54 -l3h 


(A. +rx . 

_i I_ 

30h 


36 3h -36 3h 

2 ,2 


4h 


sym 


-3h -h 
36 -3h 


4h^ 


4h‘ 


3ym 


I3h -3h 
156 22h 

4h^ 


T 

60 


36 6h - 36 
2 


2h 


sym 


0 

2 


[b] = 


4 

' 36 

j 

3h 

L 

- 36 

3h 

A^ (1 - x^. ) 


4h^ 

- 3h 

v,2 

- h 

30h 

sym 


36 

- 3h 


- 6h -h 
36 0 
6h2 
(A3.i) 


4h 



36 

6h 

- 36 

0 



"-30 

6h 30 

-6h ~ 

^3 


2h^ 

- 6h 

-h^ 

{A3 - A^) 


-6h 

0 6h 

-h" 

60 

sym 


36 

0 

60 


-30 

-6h 30 

6h 





6h^ __ 



6h 

h -6h 

0 

J 


12 

6h 

-12 

6h ; 


156 22h 54 - l3h 


aQ 

.. 2 

4h 

-6h 

2h^ 

f D h 


4h 

^ I3h -3h 

2 


3 

2h 



12 

.-6h 

” 2 X 420 . 

i 



156 22h 




sym 



4h^ ^ 

1 

sym 

4h^ 




(A3. 2) 
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12 

6h 

- 12 

6h“ 


36 3h 

- 36 3h ! 

[e] =• 

h 


4h^ 

- 6h 

2h^ 

(A +rx . ) 

- 1 

30h 

4h 

2 

- 3h -h 



12 

-6h 

sym 

36 -3h 




sym 



4h^ 

; 


4h^ 


420 


156 22h 

4h 


sym 


- I3h 


36 

6h 

-36 

0 

CM 

1 

7 


2h^ 

-6h 

-h^ 

22h 

~ 60 



36 

, 0 

4h^ 


sym 



^v2 

6h 


(A3. 3) 


[g] 


3 

60 


A3 (1 


30 h 


36 6h 


2h' 


36 


3h 

- 36 

3h 

4h2 

- 3h 

-h^ 


36 

-3h 


sym 


4h' 


sym 


36 

i 

6 


1 

-30 

6h 

30 

-6h , 

6h 

2 

-h 


-6h 

0 

6h 

-h^ 

36 

0 ; 

60 

-30 

-6h 

30 

6h 


6h^ 

-r r" 

6h 

CM 

-6h 

0 I 



12 5h 

- 12 

6h 

gQ 

4h 

- 6h 

2h^ 

2h-^ 1 


12 

-6h 


sym 


4h2 



156 22h 54 - I3h’" 

fQ h 

4h^ I3h - 3h^ 

2 420 

156 22h 


2 


sym 4h 

i J 


(A3. 4) 
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APPEITDIX - 4 


BOUNDARY CONDITIONS FOR HARMONICALLY PERTURBED FLOW 


A4-1 Simple Supports ; 

For pipe simply supported at extreme ends, the 
boundary conditions will be 

Y = 0 2 = 0 




(A4.1 ) 


For a intermediate simple support, the boundary 
conditions will be 


Y = 0 ■ 2 = 0 


(A4.2) 


A4.2 Fixed Supports : 

For fixed ends of pipe, the boundary conditions 


will be 


Y = 0 


2 = 0 


dY 

dT " 


0 


d2 

d? 


O 


A4.3 Elastic Supports, Displacement Springs 


(A4.3) 


For a displacement spring 

3 3 

^-1 + a, Y = 0, ^ +a, 2 = 0 at 5=0 
dr d 5 

^ - a . Y = O, ^ - a, 2 = 0 at5=h 

where ^ t 3 

Kd L 


(A4.4 ) 

(A4.5) 
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A4.4 Elastic Supports/ Torsional Springs : 

For pipe supported by torsional spring, 


2 

d Y 

1 

P 

ft 

II 

o 

d^Z 

I?" 


at 

o 

II 

d 

- a ^ - 0 
• ®t d5 " ° ^ 

^ + 

2 

d?^ 

Yil = ° 

at 

5 = h 


... (A4.6) 

where 

« = (A4.7) 

t El 



